Spherically symmetric solutions in f(R) gravity via 
Noether Symmetry Approach 



S. Capozziello 1 , A. Stabile 2 , A. Troisi 1 

1 Dipartimento di Scienze Fisiche, Universita di Napoli "Federico IF, and INFN, Sez. 
di Napoli, Compl. Univ. di Monte S. Angelo, Edificio G, Via Cinthia, 1-80126 - 
Napoli, Italy 

2 Dipartimento di Fisica "E.R. Caianiello" and INFN, Sez. di Napoli, Universita di 
Salerno, via S. Allende, 1-84081 - Baronissi (Sa), Italy 

Abstract. We search for spherically symmetric solutions of f(R) theories of gravity 
via the Noether Symmetry Approach. A general formalism in the metric framework 
is developed considering a point-like f(R) - Lagrangian where spherical symmetry is 
required. Examples of exact solutions are given. 

PACS numbers: 98.80.-k, 95.35,+x, 95.35.+d, 04.50.+h 
1. Introduction 

Extended Theories of Gravity have become a sort of paradigm in modern physics since 
they seem to solve several shortcomings of standard General Relativity (GR) related 
to cosmology, astrophysics and quantum field theory. The idea to extend Einstein's 
theory of gravitation is fruitful and economic with respect to several attempts which 
try to solve problems by adding new and, most of times, unjustified new ingredients 
in order to give a self-consistent picture of dynamics. The today observed accelerated 
expansion of Hubble flow and the missing matter of astrophysical large scale structures, 
are primarily enclosed in these considerations. Both the issues could be solved changing 
the gravitational sector, i.e. the l.h.s. of field equations. The philosophy is alternative 
to add new cosmic fluids (new components in the r.h.s. of field equations) which 
should give rise to clustered structures (dark matter) or to accelerated dynamics (dark 
energy) thanks to exotic equations of state. In particular, relaxing the hypothesis that 
gravitational Lagrangian has to be a linear function of the Ricci curvature scalar R, like 
in the Hilbert-Einstein formulation, one can take into account, as a minimal extension, 
an effective action where the gravitational Lagrangian is a generic f(R) function. As 
further request, one can ask for f(R) being analytical in order to recover, at least locally, 
the positive results of GR. Several studies in this sense show that cosmic dynamics at 
early [H [2] and late [3] epochs can be successfully reproduced. On the other hand, flat 
rotation curves of spiral galaxies can be fitted adopting the low energy limit of power 
law f(R) and without considering huge amounts of dark matter in galactic haloes [HE]. 
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Interesting indications have been achieved also for other f(R) functions [5] and for 
large scale structure (this issue is still matter of debate [7J, the CMBR - spectrum turns 
out to be very slightly affected if the Lagrangian shows a small deviation from the 
standard Hilbert-Einstein form [8]). Despite of these positive results, the problem is 
still open since the degeneration of viable Lagrangians has not been removed yet [9] 
and a final theory embracing the phenomenology at local and large scales or at early 
and late epochs, considering only the f(R) approach, is not available up to now (see 
[TO] [TT] for a recent discussion). On the other hand, Solar System experiments are not 
giving univocal constraints on the Parametrized Post-Newtonian limit of such theories: 
some authors claim for the fact that there is room for Extended Theories of Gravity 
considering experimental data [T2l fT3"| H%] , others claim for ruling out such theories 
with respect to GR [l~5l [T6] . The debate essentially lies on the physical meaning of the 
conformal transformations. For a recent and illuminating discussion on the argument 
see [UJ and references therein. A part the controversies, and the open discussions, 
several efforts have been done to give self-consistent formulations of f(R) gravity [IE] 
and several approaches have been pursued to find out solutions of the field equations 
coming out from such theories, both in metric and in Palatini formalism. 

In a recent paper [T9] . spherically symmetric solutions for f(R) gravity in vacuum 
have been found considering relations among functions defining the spherical metric or 
imposing a constant Ricci curvature scalar. The authors have been able to reconstruct, 
at the end, the form of some f(R) theories, discussing their physical relevance. In [20J, 
the same authors have discussed static spherically symmetric perfect fluid solutions for 
f(R) gravity in metric formalism. They showed that a given matter distribution is not 
capable of determining the functional form of f(R)- 

In this paper, we want to seek for a general method to find out spherically symmetric 
solutions in f(R) gravity and, eventually, in generic extended theories of gravity. 
Asking for a certain symmetry of the metric, we would like to investigate if such a 
symmetry holds for a generic theory of gravity. In particular for the f(R) theories. 
Specifically, we want to apply the Noether Symmetry Approach [22] in order to search 
for spherically symmetric solutions in generic f(R) theories of gravity. This means that 
we consider the spherical symmetry for the metric as a Noether symmetry and search 
for f(R) Lagrangians compatible with it. The method can give several hints toward the 
formulation of Birkhoff 's theorem (see [23] for a general formulation) for these theories 
since, up to now, there are controversial results in this direction^. 

The layout of the paper is the following. In Sec. 2, we derive and discuss the 
field equations for f(R) gravity. Sec. 3 is devoted to the construction of the point- 
like Lagrangian for a generic f(R) theory. We point out that imposing the spherical 
symmetry in the action, and then deriving the Euler-Lagrange equations, is equivalent 
to derive first the field equations and then to impose the spherical symmetry. This 

| Some authors state that the theorem is not valid in general [24] while others claims for its validity 
for specific classes of f(R) [25l[26lEZ]. 
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procedural allows to construct a Noether vector. The Noether Symmetry Approach 
to reduce dynamics is described in Sec. 4. The goal of the method is to construct a 
vector field which, contracted with the point-like Lagrangian, allows to find out, if 
existing, the conserved quantities of dynamics. Then it is possible to recast the original 
Lagrangian in a new set of variables where cyclic ones explicitly appear. The number of 
cyclic variables is equal to the number of conserved quantities. This technique reduces 
the order of derivation of the equations and simplifies the process to achieve exact 
solutions. We carry out the method for the point-like f(R) Lagrangian in spherical 
symmetry and find out some exact solutions in Sec. 5. Discussion and conclusions are 
drawn in Sec. 6. Appendix A is devoted to the field equations in spherical symmetry for 
a generic f(R)- In Appendix B, we write explicitly the PDE system, derived from the 
contraction LxC = 0, which is the existence condition for the Noether Symmetry. 



2. The f(R) gravity action and the field equations 

The action 

A= I d 4 xV=g~[f(R) + C m } } (1) 



describes a theory of gravity where f(R) is a generic function of scalar curvature, g is 
the determinant of the metric tensor and C m is the standard fluid matter minimally 
coupled with gravity. We are assuming physical units 8tcG = 1. The field equations in 
the metric approach (i.e. obtained by a variation with respect to the metric g^ u ) are 

f R (R)R^ - l -f(R)g»v - f R (R)^ + g»„nf R (R) = T™ , (2) 
df(R) 

where f R (R) = . Eqs.(T5|) are of fourth order due to the covariant derivatives of 

dR 

f R (R) and reduce to the standard Einstein ones if f(R) = R- T™ is the matter fluid 
stress-energy tensor. Defining a curvature stress - energy tensor 

T™ rv = j^r if(R) - Rfn(R)} + f R {R)^{g a ,g P u - g^g*)} , (3) 

Eqs.(T5]) can be recast in the Einstein - like form: 

1 

T 

where matter non - minimally couples to geometry through the term l/f R (R). In this 
paper, we shall seek for exact solutions in vacuum so that we can assume T™ = 0. 

By contracting with respect to the metric tensor, one can reveal the analogy of 
f R (R) with a scalar field being the trace equation 

3Df R (R) + Rf R (R) - 2f(R) = . (5) 

§ It is straightforward to show that this method works also for Friedmann- Robertson- Walker (FRW) 
metric or generic Bianchi's metrics and it is extremely useful to find out cosmological solutions 



= R^ v - -g^R = T™ v + T™/ f R (R) (4) 
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It can be seen as a Klein-Gordon equation for an effective scalar field if the identifications 

- (m dV(ip) Rf R (R) - 2f(R) 

are considered [T|. 

3. The point-like f(R) Lagrangian in spherical symmetry 

As hinted in the introduction, the aim of this paper is to work out an approach to 
obtain spherically symmetric solutions in fourth order gravity by means of Noether 
Symmetries. In order to develop this approach (the method will be outlined in Sec. 4), 
we need to deduce a point-like Lagrangian from the action ([1]). Such a Lagrangian 
can be obtained by imposing the spherical symmetry in the field action (CQ). As a 
consequence, the infinite number of degrees of freedom of the original field theory will 
be reduced to a finite number. The technique is based on the choice of a suitable 
Lagrange multiplier defined by assuming the Ricci scalar, argument of the function 
f(R) in spherical symmetry. Elsewhere, this approach has been successfully used for 
the FRW metric with the purpose to find out cosmological solutions (22J EH1 l3l] . 
In general, a spherically symmetric spacetime can be described assuming the metric : 

ds 2 = A(r)dt 2 - B(r)dr 2 - M(r)dVL , (6) 

where dQ = d6 2 + sm9 2 d(p 2 is the angular element. Obviously the conditions M(r) = r 2 
and B(r) = A~ l {r) are requested to obtain the standard Schwarzschild case of GR. This 
condition is necessary if one wants to recover the standard measure of a circumference 
when r is the radius of a circle. Our goal is to reduce the field action (JTJ) to a form with 
a finite degrees of freedom, that is the canonical action 

A= J drC{A,A',B,B',M,M',R,R') (7) 

where the Ricci scalar R and the potentials A, B, M are the set of independent variables 
defining the configuration space. Prime indicates the derivative with respect to the radial 
coordinate r. In order to achieve the point-like Lagrangian in this set of coordinates, 
we write 

A = J d 4 x^g f(R) - \{R - R) , (8) 

where A is the Lagrangian multiplier and R is the Ricci scalar expressed in terms of the 
metric 

- A" M" A'M' A 12 M' 2 A'B' B'M' 2 . , 
R = l 9 I (Q) 

AB BM ABM 2A 2 B 2BM 2 2AB 2 B 2 M M 
which can be recast in the more compact form 

A" M" 

# = + — + 2—, (10) 
AB BM ' V ; 
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where R* collects first order derivative terms. The Lagrange multiplier A is obtained 
by varying the action (jSJ) with respect to R. One gets A = f R (R). By expressing the 
determinant g and R in terms of A, B and M, we have, from Eq.(JH]), 



A = J drA 1/2 B 1/2 M^f -f R (R- R* 



A" 



M" m 



dr\A l ' 2 B l l 2 M 



f-f R (R-R*)j - ( 



AB BMJ 
f R M 



A 1 ' 2 

A'-2[^f R )M' 



A i/2 B i/2j Kb 1 / 2 

The two lines differs for a divergence term which we discard integrating by parts 
Therefore, the point-like Lagrangian becomes : 

fh' 



Al/2fR -M> 2 



A'M> - -^k^A'R' + 



2A^ 2 f RR 



B 1 / 2 



2MB 1 / 2 A 1 / 2 B 1 / 2 A 1 / 2 B 1 / 2 

R'M' - A 1/2 B 1/2 [{2 + MR)f R - Mf] 



(12) 



which is canonical since only the configuration variables and their first order derivatives 
with respect to r are present. Eq. f[T2l can be recast in more compact form introducing 
the matrix formalism : 

q'tfq' + V (13) 



C 



where q = (A, B, M, R) and q' = (A', B', M', R') are the generalized positions and 

velocities associated to C. The index "f indicates the transposed column vector. The 

d 2 C 

kinetic tensor is given by T i? - = - , - , . V = V(q) is the potential depending only on 

dq^dq'j 

the configuration variables. The Euler - Lagrange equations read 



d „ „ „ „ d . 
-V q ,C-V q C=2-[Tq> 



V q V 



= 2T g" + 2(4- V q t)i - V q V - f {V q f = (14) 

which furnish the equations of motion in term of A, B, M and R, respectively. The field 

equation for R corresponds to the constraint among the configuration coordinates. It is 

d 2 C 

worth noting that the Hessian determinant of (|T2|) . , is zero. This result clearly 



dq'idq'j 



depends on the absence of the generalized velocity B' into the point - like Lagrangian. As 
matter of fact, using a point-like Lagrangian approach implies that the metric variable B 
does not contributes to dynamics, but the equation of motion for B has to be considered 
as a further constraint equation. 

Beside the Euler- Lagrange equations (fi~4l) . one has to take into account the energy : 
E c = i- V q ,£ - C (15) 

which can be easily recognized to be coincident with the Euler-Lagrangian equation for 
the component B of the generalized position q. Then the Lagrangian f|T2l has three 
degrees of freedom and not four, as we would expected "a priori". 
Now, since the motion equation describing the evolution of the metric potential B does 
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not depends on its derivative, it can be explicitly solved in term of B as a function of 
other coordinates : 

2M 2 f RR A'R' + 2Mf R A'M' + 4AMf RR M'R' + Af R M' 2 

2AM[(2 + MR)f R — Mf] ' [ ' 

By inserting Eq. (1161) into the Lagrangian (TT2|) . we obtain a non- vanishing Hessian matrix 
removing the singular dynamics. The new Lagrangian readg[jj] 

C* = L 1/2 (17) 



with 



L = ZU = l(2 + Mfl)/ "- /M1 x 
— - M 



x [2M f RR A' R' + 2MM'(f R A' + 2Af RR R') + Af R M] . (18) 

Since —— = 0, L is canonical (L is the quadratic form of generalized velocities, A', M' 
or 

and R' and then coincides with the Hamiltonian) , so that we can consider L as the new 
Lagrangian with three degrees of freedom. The crucial point of such a replacement is 
that the Hessian determinant is now non - vanishing, being: 

3AM[(2 + MR)f R -Mf} 3 f R f RR 2 . (19) 



Obviously, we are supposing that (2 + MR)f R — Mf ^ 0, otherwise the above definitions 
of B, [Eq. (fl6l) ]. and L, [Eq. lflBl ]. lose of significance, besides we assume f RR ^ to admit 
a wide class of fourth order gravity models. The case f{R) = R requires a different 
investigation. In fact, considering the GR point - like Lagrangian needs a further lowering 
of degrees of freedom of the system and the previous results cannot be straightforwardly 
considered. From ( Tl2|) . we get : 

£ - = -2&n M ' 2 - A^Wn A ' M ' - 2A1 ' 2B1 ' 2 • < 20) 

whose Euler-Lagrange equations provide the standard equations of GR for Schwarzschild 
metric. It is easy to see the absence of the generalized velocity B' in Eq. (l20l) . Again, 
the Hessian determinant is zero. Nevertheless, considering, as above, the constraint (fl6l) 
for B, it is possible to obtain a Lagrangian with a non-vanishing Hessian. In particular 
one has : 

M rl A'M' 

BGR = JM + ^2A' (21) 



i /2 _ _ /M'(2MA' + AM') 



£*GR - L Gi? - V Jrf > ( 22 ) 



|| Lowering the dimension of configuration space through the substitution (|T6|) does not affect the 
dynamics, since B is a non-evolving quantity. In fact, introducing Eq. (|f 6j) directly into the dynamical 
equations given by (TT2|) . they coincide with those derived by (Ti8|) . 
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and then the Hessian determinant is 

J GR 



d 2 L ( 



dq-dq'j 



-1 



(23) 



which is nothing else but a non- vanishing sub- matrix of the f(R) Hessian matrix. 
Considering the Euler- Lagrange equations coming from (I2T1) and fl22|) . one obtains the 
vacuum solutions of GR, that is : 



A = k A - 



r + k\ 



B 



kok. 



2M 



A 



M = k 2 (r + k 1 ) 2 . 



(24) 



In particular, the standard form of Schwarzschild solution is obtained for k% = 0, k 2 



2GM 



and ki 



A formal summary of the field equations descending from the point -like Lagrangian 
and their relation with respect to the ones of the standard approach is given in Tab.l. 



Field equations approach 
I 

5 J d 4 x^f = 
I 

EM 9(V=9f) - 



H,„, = d, 



H = g^H. u = 



#oo = 
H rr = 
H ee = esc 2 QH W = 



H = A^Hqq - B- x H rr - 2M- 1 csc 2 6H imn = 



Point-like Lagrangian approach 

i 

5jdrC = 
I 

£V<,£ -V q C = 
E c = q;- v q ,c - c 

i 

d dC 8C_ n 

dr dA' 8A ~ U 

d 8C 8C_ rp n 

dr 8B' 8B ^ ^ C ~ U 

d 8C 8C n 

dr 8M' dM ~ U 

A combination of the above equations 



Table 1. The field-equations approach and the point-like Lagrangian approach differ 
since the symmetry, in our case the spherical one, can be imposed whether in the 
field equations, after standard variation with respect to the metric, or directly into 
the Lagrangian, which becomes point-like. The energy Ec corresponds to the 00- 
component of H^ u . The absence of B' in the Lagrangian implies the proportionality 
between the constraint equation for B and the energy function Ec ■ As a consequence, 
the number of independent equations is three (as the number of unknown functions). 
Finally it is obvious the correspondence between 88 component and field equation for 
M. The explicit form of field equations is given in App.B. 



4. The Noether Symmetry Approach 

In order to find out solutions for the Lagrangian (fl8|) . we can search for symmetries 
related to cyclic variables and then reduce dynamics. This approach allows, in principle, 
to select f(R) gravity models compatible with spherical symmetry. As a general remark, 
the Noether Theorem states that conserved quantities are related to the existence of 
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cyclic variables into dynamics [321 E2J Elj • Let us give a summary of the approach for 
finite dimensional dynamical systems. 

Let l{q l ,q l ) be a canonical, non-degenerate point-like Lagrangian where 



-Jr = ; detHij d = det 

OA 



d 2 l 

t^O, (25) 



dcfdq 



with Hij as above, the Hessian matrix related to I. The dot indicates derivatives with 
respect to the affine parameter A which, ordinarily, corresponds to time t. In our case, 
it is the radial coordinate r. In standard problems of analytical mechanics, / is in the 
form 

/ = T(q,q)-V(q), (26) 

where T and V are the "kinetic" and "potential energy" respectively. T is a positive 
definite quadratic form in q. The energy function associated with I is 

which is the total energy T + V. It has to be noted that E\ is, in any case, a constant 
of motion. In this formalism, we are going to consider only transformations which are 
point-transformations. Any invertible and smooth transformation of the "positions" 
Q l — induces a transformation of the "velocities" such that 

<3'(q) = ; (28) 

the matrix J = 1 1 dQ t /dq J \ | is the Jacobian of the transformation on the positions, and 
it is assumed to be nonzero. The Jacobian J of the "induced" transformation is easily 
derived and J ^ — » J ^ 0. Usually, this condition is not satisfied in the whole 
space but only in the neighbor of a point. It is local transformation. If one extends 
the transformation to the maximal submanifold such that J ^ 0, it is possible to get 
troubles for the whole manifold due to possible different topologies [34J. 

A point transformation Q l = Q l (q) can depend on one (or more than one) 
parameter. Let us assume that a point transformation depends on a parameter e, 
i.e. Q l = Q l {c[, e), and that it gives rise to a one-parameter Lie group. For infinitesimal 
values of e, the transformation is then generated by a vector field: for instance, as 
well known, d/dx represents a translation along the x axis, x(d/dy) — y(d/dx) is a 
rotation around the z axis and so on. In general, an infinitesimal point transformation 
is represented by a generic vector field on Q 

x = « t (q)|i- ( 29 ) 

The induced transformation (|28|) is then represented by 
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X c is called the "complete lift" of X [33]. A function /(q, q) is invariant under the 
transformation X c if 

w f ^| + (»l^ ' (31) 

where Ly^j is the Lie derivative of /. In particular, if L^cl = 0, X c is said to be a 
symmetry for the dynamics derived by I. 

In order to see how Noether's theorem and cyclic variables are related, let us 
consider a Lagrangian I and its Euler-Lagrange equations 

1\W~^ = - (32) 

Let us consider also the vector field (1301) . Contracting ( 1321) with the a"s gives 



Being 



(33) 



dX dqi dq 



cP—— = — (***-) - (34) 
dXdqi dX 1 dqi) \dX) dqi ' K ' 



from (I331) . we obtain 

s (35) 

The immediate consequence is the Noether Theorem^. 
If LyJ, = 0, then the function 

So = , (36) 
is a constant of motion. 

Remark. Eq. (1361) can be expressed independently of coordinates as a contraction of X 
by a Cartan one-form 

0, * |^ . (37) 

For a generic vector field Y = y l d/dx i , and one-form (3 = fiidx 1 , we have, by definition, 
iyP = y r Pi- Thus Eq. fl36|) can be written as 

i*0i = S . (38) 
By a point-transformation, the vector field X becomes 

x =(«'»^ + (IH^' (39) 

% In the following, with abuse of notation, we shall write X instead of X c , whenever no confusion is 
possible. 
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We see that X' is still the lift of a vector field defined on the "space of positions". If X 
is a symmetry and we choose a point transformation such that 

= i^l, (40) 

. (41) 

Thus Q 1 is a cyclic coordinate and the dynamics can be reduced [32], [33] . 
Remarks: 

(i) The change of coordinates defined by (140]) is not unique. Usually a clever choice is 
very important. 

(ii) In general, the solution of Eq. (l4"0~j) is not well defined on the whole space. It is local 
in the sense explained above. 

(iii) It is possible that more than one X is found, say for instance X 1; X 2 . If they 
commute, i.e. [X!,X 2 ] = 0, then it is possible to obtain two cyclic coordinates by 
solving the system 

i^dQ 1 = 1; i^dQ 2 = 1; i^dQ 1 = 0; i ^ 1; i^JQ 1 = 0; % ± 2 . (42) 

The transformed fields will be d/dQ 1 , d/dQ 2 . If they do not commute, this 
procedure is clearly not applicable, since commutation relations are preserved by 
diffeomorphisms. Let us note that X 3 = [X 1; X 2 ] is also a symmetry, indeed, being 
L^J = L^L^ 2 l — L^ 2 L^J = 0. If X 3 is independent of Xi, X 2 , we can go on 
until the vector fields close the Lie algebra. The usual way to treat this situation is 
to make a Legendre transformation, going to the Hamiltonian formalism and to a 
Lie algebra of Poisson brackets. If we look for a reduction with cyclic coordinates, 
this procedure is possible in the following way: 

• we arbitrarily choose one of the symmetries, or a linear combination of 
them, and get new coordinates as above. After the reduction, we get a new 
Lagrangian /(Q); 

• we search again for symmetries in this new space, make a new reduction and 
so on until possible; 

• if the search fails, we try again with another of the existing symmetries. 

Let us now assume that I is of the form (126]) . As X is of the form (130]) . L^l will 
be a homogeneous polynomial of second degree in the velocities plus a inhomogeneous 
term in the q 1 . Since such a polynomial has to be identically zero, each coefficient 
must be independently zero. If n is the dimension of the configuration space, we get 
{1 + n(n + l)/2} partial differential equations (PDE). The system is over determined, 
therefore, if any solution exists, it will be expressed in terms of integration constants 
instead of boundary conditions. It is also obvious that an overall constant factor in the 
Lie vector X is irrelevant. In other words, the Noether Symmetry Approach can be used 
to select functions which assign the models and, as we shall see below, such functions 



iy.dQ 1 = 1 ; i yi dQ i 

we get 
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(and then the models) can be physically relevant. This fact justifies the method at least 
a posteriori. 

5. The Noether Approach for f(R) gravity in spherical symmetry 

Since the above considerations, if one assumes the spherical symmetry, the role of the 
affine parameter is played by the coordinate radius r. In this case, the configuration 
space is given by Q = {A, M, R} and the tangent space by TQ = {A, A', M, M', R, R'}. 
On the other hand, according to the Noether theorem, the existence of a symmetry for 
dynamics described by the Lagrangia njfrS]) implies a constant of motion. Let us apply 
the Lie derivative to the ( fl8l) . we havdZJ: 

L X L = a- V ? L + o/- V q >L = f a - V g L + 2 f V ff a) t q, (43) 
which vanish if the functions a satisfy the following system (see App.B for details) 

a ■ V g L + 2(V,a)*L = — + 2^t im = . (44) 

OQi "Ik 

Solving the system (j44j) means to find out the functions ojj which assign the Noether 
vector. However the system (|44p implicitly depends on the form of f(R) and then, by 
solving it, we get also f(R) theories compatible with spherical symmetry. On the other 
hand, by choosing the f(R) form, we can explicitly solve ( )44l) . As an example, one finds 
that the system (T44]) is satisfied if we chose 

f(R) = f R s a = (a u a 2 ,a 3 ) = ({3 - 2s)kA, -kM, kR\ (45) 

with s a real number, k an integration constant and /o a dimensional coupling constant^]- 
This means that for any f(R) = R s exists, at least, a Noether symmetry and a related 
constant of motion So : 

S = a ■ V g 'L = 



= 2skMR 2s - 3 [2s + (s- l)MR][(s - 2)RA' - (2s 2 - 3s + l)AR'} . (46) 

A physical interpretation of So is possible if one gives an interpretation of this quantity 
in GR. In such a case, with s = 1, the above procedure has to be applied to the 
Lagrangian (1221) . We obtain the solution 

a GR = (-kA, kM) . (47) 

The functions A and M give the Schwarzschild solution (|24[) . and then the constant of 
motion acquires the form 
2GM 



c 2 



( 4 * 



+ From now on, q indicates the vector (A, M, R). 

* The dimensions are given by i? 1_s in term of the Ricci scalar. For the sake of simplicity we will put 
/o = 1 in the forthcoming discussion. 
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In other words, in the case of Einstein gravity, the Noether symmetry gives as a conserved 
quantity the Schwarzschild radius or the mass of the gravitating system. 

Another solution can be find out for R = Rq where Rq is a constant (see also [19|). 
In this case, the field equations ^ reduce to 

Rfiv + k g„u = , (49) 
where k = —-f(R )/fji(Ro). The general solution is 

A(r) = -l- = l + ^ + ^V, M = r 2 (50) 
B[r) r 12 

with the special case 

A(r) = -±- = l + ^, M = r 2 , R = 0. (51) 
B[r) r 

The solution (1501) is the well known Schwarzschild-de Sitter one which is a solution in 

most of modified gravity theories. It evades the Solar System constraints due to the 

smallness of the effective cosmo logical constant. However, other spherically symmetric 

solutions, different from this, are more significant for Solar System tests. 

In the general case f(R) = R s , the Lagrangian ( 1181) becomes 

_ sR 2s - s [2s + (s- I) MR] 
L = x 

M 

x [2(s - 1)M 2 A'R' + 2MRM'A' + 4(s - 1)AMM'R' + ARM 12 } , (52) 
and the expression ffl6|) for B is 

s[2(s - \)M 2 A!R' + 2MRM' A' + 4(s - \)AMM'R' + ARM' 2 } 
~ 2AMR[2s +(s- l)MR} ' ' 

As it can be easily checked, GR is recovered when s = 1. 

Using the constant of motion (1461) . we solve in term of A and obtain 

4s 2 -9s+5 

A = R —3 {ki + En / ; : : ; : r \ (54) 

I °J 2ks(s - 2)M[2s + (s - I) MR] > V ; 

for s / 2, with h\ an integration constant. For s — 2, one finds 

12kr 2 (4 + r 2 R)RR' ' v ; 

These relations allow to find out general solutions for the field equations giving the 
dependence of the Ricci scalar on the radial coordinate r. For example, a solution is 
found for 

s = 5/4, M = r 2 , R = 5r- 2 , (56) 
obtaining the spherically symmetric metric 

ds 2 = -^=(k 2 + hr)dt 2 - - ( i + fa ) dr 2 - r 2 dQ , (57) 
32S 

with A;2 = r- It is worth noting that such exact solution is in the range of s values 

ruled out by Solar System observations, as pointed out in 
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6. Discussion and Conclusions 

In this paper, we have discussed a general method to find out exact solutions in Extended 
Theories of Gravity when a spherically symmetric background is taken into account. In 
particular, we have searched for exact spherically symmetric solutions in f(R) gravity by 
asking for the existence of Noether symmetries. We have developed a general formalism 
and given some examples of exact solutions. The procedure consists in: i) considering 
the point-like f{R) Lagrangian where spherical symmetry has been imposed; ii) deriving 
the Euler-Lagrange equations; Hi) searching for a Noether vector field; iv ) reducing 
dynamics and then integrating the equations of motion using conserved quantities. 
Viceversa, the approach allows also to select families of f{R) models where a particular 
symmetry (in this case the spherical one) is present. As examples, we discussed power 
law models and models with constant Ricci curvature scalar. However, the above method 
can be further generalized. If a symmetry exists, the Noether Approach allows, as 
discussed in Sec. 4, transformations of variables where the cyclic ones are evident. This 
fact allows to reduce dynamics and then to get more easily exact solutions. For example, 
since we know that f{R) = R s - gravity admit a conserved quantity, a coordinate 
transformation can be induced by the Noether symmetry. We ask for the coordinate 
transformation : 



L = L(q,q) = L(A, M, R, A', M' , R') -»• L = L(M, R, A', M',R') , (58) 

for the Lagrangian (fl8|) , where the Noether symmetry, and then the conserved quantity, 
corresponds to the cyclic variable A. If more than one symmetry exists, one can find 
more than one cyclic variables. In our case, if three Noether symmetries exist, we 
can transform the Lagrangian L in a Lagrangian with three cyclic coordinates, that is 
A = A(q), M = M(q) and R = R{q) which are function of the old ones. These new 
functions have to satisfy the following system 

(3-2s)A§-M§ + R§ = l, 

(59) 

(3-2s)A& -M^ + R 9 ^ = 0, 
with i = 2, 3 (we have put k — 1). A solution of (I591) is given by the set (for s ^ 3/2) 
A = J^U + F A (A^M^ A , A— M^ A ) 

(60) 

ft = F^A^M 1 *, A^M 6 ) 
and if s = 3/2 

' A = -In M + F A (A)G A (MR) 

(61) 

ft = Fi(A)Gi(MR) 

where F A , F i: G A and Gi are arbitrary functions and r] A , rji, £ A and & integration 
constants. 
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These considerations show that the Noether Symmetries Approach can be applied 
to large classes of gravity theories. Up to now the Noether symmetries Approach has 
been worked out in the case of FRW- metric. In this paper, we have concentrated 
our attention to the development of the general formalism in the case of spherically 
symmetric spacetimes. Therefore the fact that, even in the case of a spherical symmetry, 
it is possible to achieve exact solutions seems to suggest that this technique can represent 
a paradigmatic approach to work out exact solutions in any theory of gravity. At this 
stage, the systematic search for exact solution is well beyond the aim of this paper. A 
more comprehensive analysis in this sense will be the argument of forthcoming studies. 
A final comment deserves the possible relevance of this approach for the above mentioned 
Birkhoff- Jensen theorem. The validity of such a theorem is crucial in every theory of 
gravity, due to the fact that it is directly related to the physical properties of self- 
gravitating systems (stability, stationarity, etc.). The results presented in this paper 
point out that it does not hold in general for the specific f(R) theories considered. 
However, the above technique could be a good approach to select suitable classes of 
theories where such a theorem holds. 

Appendix A. The f(R) field equations in spherical symmetry 

The field equations (Tab.l) in spherical symmetry, derived from the variational principle 
of the action (CQ), are 



#00 = 2A 2 B 2 Mf + {BMA' 2 - A[2BA'M' + M{2BA" - A'B')]}f R + 



+ (-2A 2 MB'R' + AA 2 BM'R' + AA 2 BMR")f RR + 



+AA 2 BMR' 2 f RRR = 0, 



(A.l) 



H„ = 2A 2 B 2 M 2 f + (BM 2 A' 2 + AM 2 A'B' + 2A 2 MB'M' + 2A 2 BM' 2 + 



2 ABM 2 A' - AA 2 BMM")f R + (2ABM 2 A'R' + 



+AA 2 BMM'R')f RR = 



(A.2) 



Hob = 2AB 2 Mf + (AAB 2 - BAM' + AB'M' - 2ABM")f R + 



+ (2BMA'R' - 2AM B'R' + 2ABM'R' + AABMR")f RR + 
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+4ABMR' 2 f RRR = 0, (A.3) 



H w = sin 2 6Hee = 



(A.4) 



The trace equation is 



H = g^H^ = AAB 2 Mf - 2AB 2 MRf R + 3(BMA'R' - AMB'R' + 



+2ABM'R' + 2ABMR")f RR + 6ABMR' 2 f RRR = 



(A.5) 



Appendix B. The Noether vector 



The system (|44)) comes out from the condition LxL = for the existence of the Noether 
symmetry. Considering the configuration space q — (A , M , R ) and defining the Noether 
vector components a = («i , a 2 , a 3 ), the system (H4|) assumes the explicit form 



da 3 



(B.l) 



M 



2 + MR)a 3 f RR -'^f R ]f R + 

2Ada- 
+ ' 



+ HvM + dM 



M dM J \M dM J 



(B.2) 



e(^ + 2A^)/ JUl = 



«9i? 



«9i? 



(B.3) 



a 3 + M^ + 2A^)f RR + 



dM 



dA 



+ 



<9a 2 da i A <9a 2 
+ ^4 + ~M~dA 



f R ] 



(B.4) 



[:V/(2 + :V/./?)n,./ /l », v -2no./ /,]./»/,. -K + 



/ (9«i Cao oa 3 \ i 

+ ^2a 2 + + 2A-^- + M-^ja 3 f RR + M f RRR j = (B.5) 



da-. 



da 



Spherically symmetric solutions in f(R) gravity via Noether Symmetry Approach 16 



Y<9«i t A On 2 

da 3 r dai da: 



2A[(2 + MR)a 3 f RR - (f - Rf R )a 2 ]f RR + + --^ )f R 



with the condition £ = (2 + MR)f R — Mf ^ 0, otherwise the Hessian of Lagrangian 
L (Fl8|) is vanishing. 
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